module real_space_diagnostic
use mod_2dflu
implicit none
public
contains
!!--------------------
subroutine initialize_rspace_diag()
!
endsubroutine initialize_rspace_diag()
!!
!!
subroutine force_inj
  use mod_2dflu
  implicit none
  integer :: i1,i2
  double precision :: rk2,x,y
!!!
  fw = 0.0d0; fu = 0.0d0
!
  do i2 = 1,n2
    do i1 = 1,n1
      x = i1*dx
      fu = fu + famp*uky(i1,i2)*dsin(kf*x)
      fw = fw + famp*kf*omega(i1,i2)*dcos(kf*x)
    enddo
  enddo
  fw  = fw /(dble(n1)*dble(n2))
  fu  = fu /(dble(n1)*dble(n2))
  !! --------------------- 
end subroutine force_inj
!!-------------------------
subroutine calc_stfunc()
endsubroutine calc_stfunc
!-------------------------
end module real_space_diagnostic
